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Abstract 

We discuss pairing correlations in weakly bound neutron rich nuclei, by using 
the coordinate-space Hartree-Fock-Bogolyubov approach which allows to take 
properly into account the coupling to particle continuum. We show that the 
additional pairing binding energy acts against a development of an infinite 
rms radius, even in situations when an i=0 single-particle orbital becomes 
unbound. 
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An essential element in understanding the physics of weakly bound fermion systems, 
is the residual coupling between bound states and continuum. With decreasing particle 
separation energy, the zone of correlated states around the Fermi surface must incorporate 
the phase space of particle continuum. One of the most interesting discoveries in loosely 
bound nuclei is the neutron halo In this low density, high-isospin state of matter, the 
two-particle correlations, like the pairing modified by the effects of virtual scattering to the 
continuum, are expected to play a crucial role. 

Nuclear pairing correlations are important throughout the periodic table, providing an 
extra binding for ordinary even-even and odd nuclei, but they become essential for loosely 
bound nuclei, because without them the systems in question cease to be bound. For their 
description, two main elements have to be combined. Firstly, the residual interaction which is 
expected to be different in the nuclear exterior and in the nuclear interior, should include the 
density dependence in the self-consistent framework. Secondly, the residual coupling between 
bound and continuum states must be treated consistently. Both elements are present in the 
Hartree-Fock-Bogolyubov (HFB) coordinate-space, mean-field approach using the Skyrme 
energy density functional 0. 

The asymptotic behavior of halo nuclei, was discussed both for 2-body and 3-body 
systems For 2-body systems, the root-mean-square radius (rms) of a low-angular- 

momentum (£=0,1) neutron orbital in an average potential diverges when its energy ap- 
proaches the continuum. Rigorous results for 3-body systems have been obtained neglecting 
the antisymmetrization. The strong divergencies seen in 2-body systems are now seldom. 
When all 2-body subsystems of the 3-body system are unbound, the rms radius diverges at 
most logarithmically with vanishing binding energy 0. When one subsystem has a bound 
state with energy close to the total energy of the 3-body system, the system behaves as one 
particle in the common field of the 2-body state Finally, when at least two subsystems 
have bound states, the ground state (g.s.) rms radius always converges [^. 

The pairing near the neutron drip line was studied first in a schematic 3-body model of 
two interacting nucleons together with a structureless core The exponential decay of 
both the single-particle (s.p.) density and the 2-particle density for valence neutrons was 
found to be governed rather by the 2-particle energy of the correlated pair than by the s.p. 
binding energy. This analysis is insufficient because just near the drip lines the particle- hole 
(p-h) and particle-particle (p-p) channels of the effective interaction are strongly coupled, 
and the non-self-consistent treatment of the pairing interaction based on computing pairing 
correlations after determining the s.p. basis is unjustified. 

In light nuclei, many-body correlations may essentially modify the mean-field picture 
|^,TH|. For that, different cluster models have been applied [rT,T2|, and the 3-body dynam- 



ics has been investigated 11,13-15]. Much progress has been made in the treatment of the 



Pauli exclusion principle ||16|, which is needed to remove spurious components of the halo 



wave function even at small densities. The whole problem is subtle because the proper anti- 
symmetrization of the wave function easily spoils many advantages of the cluster description. 
It was found that the 3-body asymptotics of the cluster wave function is given by the expo- 



nential fall-off with a decay constant /t=y 2m£'(3)/^, where i?(3) is the binding energy of the 
system relative to the 3-body break-up threshold. If a bound 2-body subsystem is present, 
the asymptotics is given by a binary channel asymptotics [|l^. 

In heavier nuclei, the concept of self-consistent mean-field is better justified. Moreover, 
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the one-body densities are always governed by the so-called natural orbitals (the eigenstates 
of the exact one-body density matrix), and the mean- field approximation can be understood 
as a method to analyze properties of the natural orbitals. Strong pairing correlations do not 
allow for the decoupling of weakly bound neutrons from the rest of the system [0,0 . It is 



therefore a crucial question whether the neutron halo phenomenon can happen close to the 
neutron drip-line in heavy nuclei. 

In the case of no pairing, irrespective of whether the number of neutrons is even or odd, 
the g.s. neutron mean-field asymptotic density equals p(r)~exp(— 2Kr)/r^ for k=\/ —2me/h 
[e is the s.p. energy of the least bound £=0 neutron), and has a long tail with the rms radius 
Tj-^s^l/ K^oo for e^O. 

In the presence of pairing, we have to separately consider the even-A^ and odd-A^ isotopes. 
For N even, the asymptotics of the HFB neutron density is given by the binding energy 
relative to the 2-body break-up threshold E(2)- Therefore, the decay constant becomes 



K=^2mE{2)/h for i?(2)=min[£'^] — A, where min[ii^^] is the lowest discrete {E^<—X) 
quasi-particle (q.p.) energy, and A is the chemical potential. The HFB q.p. energies are 



fairly close to the canonical q.p. energies pi 



E, = y'(e- - A)2 + (A-)^, (1) 

where e^^"^ and A^*^"^ are the diagonal matrix elements of the p-h and p-p mean fields in the 
canonical basis. Hence, a canonical state at the Fermi surface (e^'^^^A) gives the asymptotic 
density for £'(2)=A^'^°— e^*^"^, which is larger by A^*^"^ than that in the no-pairing approxima- 
tion. Therefore, in even-iV nuclei, the paired densities decrease faster than the unpaired 
densities, which is the effect called here the 'pairing anti-halo effect'. 



In odd-A^ nuclei, the neutron densities are given by one-q.p. states p2[, which gives the 
asymptotic decay corresponding to E(^2) = — iiiin[i?^] — A. This threshold energy goes to zero 
when the minus Fermi energy approaches the pairing gap of a state at the Fermi surface. 
Therefore, the pairing correlations do not prevent the neutron radii of odd-A^ nuclei from 
diverging in the limit of A~ gcan_^_^can_ 

Below we give a numerical example illustrating the mechanism of the pairing anti-halo 
effect in even-A^ nuclei, which takes the most dramatic form for e^^'^^A-^O, because then the 
no-pairing tail of the neutron density gives a divergent rms radius, while that for the paired 
density stays finite. Before that we must, however, consider that in weakly bound nuclei the 
situation is more complicated, because for A— s>0 there can be no discrete HFB quasiparticle 
state at all, and we have to take into account the HFB continuum explicitly. To this end we 
express the asymptotic density, coming from the q.p. continuum, in the form of the integral: 

p(r)~/ dE^^expi-2KiE)r), (2) 



where K{E)=^2m{E — X)/h, and A[E) is the amplitude of the exponential tail correspond- 
ing to the quasiparticle having the q.p. energy E. It is obvious that at r^oo the upper 
limit of the integral is not important, and the asymptotic behavior of the density is 
determined by amplitudes A{E) near the threshold of E=—X. 

These amplitudes are determined by the pairing interaction and cannot be inferred from 
any simple geometric considerations. We can expect that A{E)—*0 for k,{E)^0, because 
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in this limit, the lower HFB q.p. wave functions are completely delocalized |T8[, and thus 
become decoupled from the pairing field inside the nucleus. It is instructive to realize that 
a power-law dependence of A)" gives for A=0 the density vanishing according 

to a power law too: p(r)~r^'^^^". At any small non-zero value of A, the latter power law 
holds up to r~l/K(— A), but then turns over to the standard exponential dependence of 
exp(— 2k(— A)r). 

The two-neutron drip line (A~0) in heavier nuclei will not probably be reached before 
soon, however, the experimental data for neutron halos exist in light nuclei [^. Fig. [^(a) 
compares the calculated one-neutron separation energies, obtained using the spherical HFB 
code of Ref. with the experimental and systematics values in heavy carbon isotopes 



The Skyrme SLy4 interaction in the p-h channel together with the zero-range, 



density-dependent (surface) pairing p-p interaction, as described in Refs. [pT| , p6[ , has been 



used in these calculations. The overall agreement is encouraging, if one considers that the 
spherical symmetry is assumed and, moreover, correlations beyond the mean-field are very 
important in these light systems. 

The HFB code has been improved by employing the analytical wave functions in the 
asymptotic region, which allowed for solving the HFB equations in a very large box, and 
thus obtaining an arbitrarily dense discretized-continuum spectrum. The box of -Rbox=125 fm 
has been used in the present study, which gives the lowest discretized-continuum HF state at 
e~(7r^/;,^)/(2mi?^Qj^)~0.013 MeV. Therefore, we may safely treat the mean-field continuum 
phenomena for the Fermi energies —A at least several times larger than this limiting value 
of 0.013 MeV. We have also performed test calculations in the box of i?box=400fm, for 
which this limiting value is an order of magnitude smaller, and we found no modifications 
of conclusions reached in the present study. 

In Fig. |I|(b) the spherical self-consistent Hartree-Fock (HF) neutron s.p. energies are 
shown for Z=Q as functions of the neutron number A^. The carbon isotopes are special in 
that the least bound spherical neutron orbital at the neutron drip line has £=0, i.e., infinite 
rms radius should appear in the no-pairing scenario, had this orbital approached the zero 
binding energy. Similarly, an infinite rms radius would be given by an unpaired weakly 
bound c?e/ormed orbital, had this orbital contained even a small 2si/2 component. 

In self-consistent calculations within a given effective interaction, it is extremely difficult 
to change the relative positions of the s.p. orbitals without spoiling other global character- 
istics of the weakly bound nuclei. Therefore, in order to study dependence of the results 
on the 2si/2 binding energy, we have performed the HFB calculations for fixed (i.e., non- 
self-consistent) s.p. wave functions of the Poschl-Teller-Ginocchio (PTG) potential (see the 
description and references given in Ref. |TE]). Parameters of the PTG potential have been 
fixed at A=5, s=0.085fm~^. The depth parameter of the PTG potential for the si/2, P1/2, 
P3/2, (i3/2, and ^5/2 channels is z/£j=3.251, 2.814, 3.126, 2.747, and 3.176, respectively. At 
A^=16, these values reproduce the HFB+SLy4 neutron spectrum of Fig. ^b) exactly. Using 
the PTG potential we may now move the 2si/2 orbital up or down, by changing the value 
of the depth parameter z/si/2. 

Figure |(a) shows the s.p. energy of the 2si/2 orbital as function of the potential depth 
With decreasing the 2si/2 orbital approaches the zero binding energy at z/si/2=3 
(full circles), then turns into a virtual state (dashed line), and finally reappears at 1/^1/2=2.857 
as a very broad resonance (open circles). The rms radius of this orbital becomes infinite 
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at e=0, and hence the neutron PTG rms radius shown in Fig. |5(b) for A^=16 (no pairing) 
becomes infinite too. On the other hand, the HFB rms radius which is also shown in 
Fig. 0(b), is completely unaffected by the transition of the 2si/2 level through the point of 
no s.p. binding. This constitutes a vivid illustration of the pairing anti-halo effect. Detailed 
behavior of the rms radius can be traced back to the behavior of the 2si/2 q.p. energies (|1]) 
induced by the canonical s.p. energies e^^"^ and pairing gaps A^'^'^ shown in Figs. H(a) and 
(c), respectively. 

The HFB calculations have to be performed not at a fixed particle number of A^=16, 
but at a fixed Fermi energy of A=— 0.056 MeV. This value of A corresponds to A^=16 at 
z/si/2=3.251. At smaller values of 2/^1/2, the 2si/2 orbital depopulates, and the particle 
number cannot be kept constant for any negative value of A. Values of the average neutron 
number {N) at A=— 0.056 MeV (see Fig. 0(c)), change from 16 to about 14.7 in the studied 
interval of i'si/2- Nevertheless, the canonical-basis occupation factor v"^^^ of the 2si/2 
canonical state (open squares in Fig. 0(c)), is far from zero and decreases only to about 0.16 
at the smallest value of z/si/2=2.778. We have also performed calculations for fixed A^=14 
(this value can be kept constant with the varying position of the 2si/2 orbital), and we 
obtained exactly the same qualitative results as for fixed A and a variable neutron number. 
The occupation of the 2si/2 canonical state for A~0 is entirely due to virtual pair excitations. 

In order to visualize the origins of the pairing anti-halo effect, in Fig. ^ are shown the 
contributions, N{E) and r^{E), of the i=0 continuum HFB q.p. states to the particle number 
and radius squared, respectively. They are defined as the corresponding space integrals of 
the lower components of the q.p. wave functions squared, |v?2(-E', r)p. When the 2si/2 
orbital is well bound (see panels (a) and (d)), the largest contributions to the particle number 
(N) = dE N (E) , and to the radius squared {r'^)=J^^dEr'^{E), come from narrow peaks 
located near the corresponding s.p. energies. However, when the 2si/2 orbital is pushed 
across 6=0, and up into the continuum, the peaks in N{E) and r'^{E) do not at all approach 
£^=0, but stay at a fairly large values of ii^~A^'^°~2 MeV. In this sense, the large pairing gaps 
prevent the largest contributions to the rms radius from entering the region of divergencies 
at ii^~0, where the lower components of the HFB q.p. wave functions may develop long tails. 

The sequence of curves presented in Figs. 0(d)-(f) is particularly instructive, because 
it shows that the 2si/2 orbital gives either large or small contribution to the rms radius, 
depending on whether it is bound or resonant, but this contribution is always located around 
ii^~2MeV and simply disappears with 2si/2 being pushed up. On the other hand, the 
low-energy continuum contributions, described by the asymptotic density matrices given in 
Eq. (0), are almost independent of the 2s 1/2 position, and give well- pronounced, low-energy 
contributions to the rms radii. However, at any small but non-zero value of A, the curve 
r'^{E) bends over at small E, and goes to zero, thus giving a finite, non-diverging values of 
the HFB rms radii. 

In conclusion, we have shown that pairing correlations significantly modify the asymp- 
totic properties of one-body neutron densities p(r) in weakly bound neutron rich nuclei. The 
neutron rms radii in even-iV nuclei are predominantly determined by two types of quasi- 
particles: (i) those near the threshold of E=—X, and (ii) those corresponding to the i=0 
s.p. orbitals. We have checked by numerical calculations that: (i) the density of the former 
quasiparticle states vanishes at E=—\, which is due to a delocalization of the quasiparticle 
wave functions, and (ii) the latter quasiparticles have energies E larger than the correspond- 
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ing pairing gaps, and thus cannot give diverging rms radii. As a result, and contrary to the 
no-pairing approximation, the rms neutron radii in even-N nuclei never diverge, even in the 
limit of A=0. 

In this work, we have analyzed the HFB asymptotics corresponding to the one- 
neutron densities. The 3-body asymptotics of two-neutron densities cannot be consis- 
tently described in the HFB. One should notice, however, that the break-up threshold 
for two-neutron emission is expected to be higher than the two-neutron separation energy 
S2N=E^''^'"^\N-2, Z)-E(^''^''\N, Z) by the amount of the 2-particle energy of the corre- 
lated neutron pair. Hence, one expects that the pairing correlations will also act against 
the two-neutron halo effect in heavy even-N systems, by preventing the neutron radii from 
rapid growth in the limit of S2n^0. 
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FIGURES 




FIG. 1. One-neutron separation energies Sn (a) and the spherical s.p. energies (b) in heavy 
carbon isotopes. Spherical self-consistent HFB and HF calculations have been performed for the 
Skyrme SLy4 interaction |^5|. Experimental (EXP) and systematics (SYS) separation energies are 
taken from Ref. \24\. 
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FIG. 2. Results of the spherical HFB calculations performed for the model PTG s.p. spectrum 
with varying positions of the 2si/2 level. Panel (a) shows the HFB canonical s.p. energies (open 
squares), the PTG s.p. energies of the bound 25^/2 states (full circles) and the real parts of the 
2si/2 resonances (open circles). The HFB neutron rms radius (b), average pairing gap (A) (c), 
and average neutron number (N) (d) (full squares in panels (b)-(d)), were calculated for a fixed 
value of the neutron Fermi energy of A=— 0.056 MeV. The corresponding canonical pairing gap (c) 
and average neutron number (d) are shown with open squares. The no-pairing PTG rms radius in 
panel (b) corresponds to N=16. 
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FIG. 3. Contributions N{E) and r^{E) to the average neutron number and radius squared, 
respectively, calculated for the HFB q.p. states, as functions of the q.p. energy E. The HFB 
calculations have been performed at a fixed Fermi energy of A=— 0.056 MeV. The curves correspond 
to calculations performed for different positions of the 2sx/2 PTG states, with the s.p. energies 
indicated explicitly [see also Fig. |2|(a)]. 
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